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ABSTRACT 

We study the contribution of a stochastic background (SB) of primordial mag- 
netic fields (PMF) on the anisotropies in temperature and polarization of the 
cosmic microwave background radiation (CMB). A SB of PMF modelled as a 
fully inhomogeneous component induces non-gaussian scalar, vector and ten- 
sor metric linear perturbations. We give the exact expressions for the Fourier 
spectra of the relevant energy-momentum components of such SB, given a 
power-law dependence parametrized by a spectral index ng for the magnetic 
field power spectrum cut at a damping scale ku. For all the values oi ub con- 
sidered here, the contribution to the CMB temperature pattern by such a SB 
is dominated by the scalar contribution and then by the vector one at higher 
multipoles. We also give an analytic estimate of the scalar contribution to the 
CMB temperature pattern. 

Key words: Cosmology: cosmic microwave background - Physical data and 
processes: magnetic fields. 



1 INTRODUCTION 

The origin of the large scale magnetic fields observed 
is an issue of grea t importance in astrophysics (see 
ISubramanianI ([2006') for a review) . Primordial magnetic 
fields (PMF) generated in the early Universe are a possi- 
ble explanation of large scale magnetic fields in clusters 
of galaxies and galaxies and might have left an imprint 
in the anisotropy pa ttern of the co smic microwave back- 
ground (CMB) (see iDurreij l|2007l ) for a review) . 
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A stochastic background (SB) is the simplest 
way to model these random PMFs in an isotropic 
and homogeneous background. A SB of PMF is 
modelled as a fully inhomogeneous component and 
its energy momentum tensor (EMT) - quadratic 
in the magnetic fields - is considered at the same 
footing as linear inhomogeneities in the other compo- 
nents and linear metric fluctuations. A SB of PMF 
generates independent modes for all kinds of linear 
perturb ations: there has b e en several studies for 
scalar and Leo (2000); Kahnias hvili and Ratral 

(|2006l ): [Ciovannini and Ku nzc (20QS) and refer- 
ences t herein ; Finelli, Paci fc Paoletti 2008 - hence- 
forth IFPPI -), vector (|Subramanian and Barrowl 
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19981: 



Seshadri and Subramaniaiil 



200C: 



2001 



Mack. Kahniashvili and Kosowskvl 20021: Lewisj 



2004) an d tensor (Durrcr. Ferreira and Kahniashvili' 



Mack. Kahniashvih and Kosowsky 



Caprini. Durrer and Kahniashvih 120041 ') perturbat 



2002; 



in presence of a SB of PMF. Limits on the amplitude 
and spectral index of the PMF were also obtained by an 
exploration of a flat ACDM model in presence of such 
SB (.Yamazaki. Ichiki. Kaiino and Mathews 2006.). 

We study the problem by a numerical Eistein- 
Boltzmann cod e exte nding the results obtained in our 
previous work (|FPP|) . The study of the impact of a 
SB of PMF on CMB anisotropics requires a detailed 
study of the initial conditions for fluctuations and of 
the power spectra of the EMT of the SB of PMF. 
Our paper improves previous results in these two as- 
pects. In this paper we obtain the Fourier spectra 
of the relevant vector and tensor energy-momentum 
compon ents o f the SB of PMF along the procedure 
used in IFPPI for the scalar components. As shown in 
IFPPI . by solving exactly the convolution integrals for 
a sequence of values for the spectral slope n_g which 
p arametrizes the PMF power spectrum, pre vious results 
of lMack. Kahniashvili and Kosowskvl l|2002l ) may be sig- 
niflcantly improved. 

With these improved correlators we then investi- 
gate the impact of a stochastic background of primordial 
magnetic fields on scalar, vector and tensor cosmolog- 
ical perturbations and in particular on CMB temper- 
ature and polarization anisotropics. Our results show 
that it is very important to study also vector perturba- 
tions since these dominate at high I over the scalar ones 
for any slope of the spectrum of PMF. 

Our paper is organized as follows. In Section 2 we 
introduce our conventions for a non-helical SB and the 
scalar, vector and tensor decomposition of its EMT: 
we give our exact results for the relevant objects for 
a set of n_B leaving the details in Appendices A, B, 
and we show how our exact results improve on pre- 
vious results. The set of ns now includ es als o values 
n_B < —3/2 which were not studied in iFPPl . Section 
3 presents the decomposition of metric perturbations 
and sections 4,5,6 present the study of scalar, vector 
and tensor cosmological perturbations in presence of a 
SB of PMF, respectively. In Section 7 we present the 
result s obtained by our modified version of the CAMB 
code (|Lewis.Challinor and Lasenblvl2000l ). We conclude 
in Section 8. 



2 STOCHASTIC BACKGROUND OF 
PRIMORDIAL MAGNETIC FIELDS 

Since the EMT of PMF at homogeneous level is zero, at 
linear order PMFs evolve like a stiff source and there- 
fore it is possible to discard the back-reaction of gravity 
onto the SB of PMF. Before the decoupling epoch the 



electric conductivity of the primordial plasma is very 
large, therefore it is possible at lowest order to con- 
sider the infinite conductivity limit, in which the in- 
duced electric field is zero. Within the infinite conduc- 
tivity limit the magnetic field amplitude scales simply 
as B(x,t) = B(x)/a(r)2 and Q the EMT of a SB of 
PMF is: 

B^(x) 



To = 



-PB = 



87ra4(T) 



= 



47r a4(r) 



S'(x)„ 



5]-B,(x)B'(x; 



(2.1) 
(2.2) 
(2.3) 



The two point correlation function in the Fourier space 
for fully inhomogeneous fields is: 



{i3.(k)B;(k')) = (2^)^5(k-k')(<5., 



PB{k) 



(2.4) 



Where Psik) is the spectrum of PMF parametrized as: 
^ (2.5) 



and fc, is a reference scale. PMF are damped on small 
scales by radiation viscosity. We model this damping 
introducing a sharp cut-off in the PMF power spectrum 
at a damping scale called kn ■ The relation between the 
amplitude of PMF power spectrum and the amplitude 
of the field itself is: 



= — / dkk'PBik) 
Solving the integral above we obtain: 



A 



27r2(ris+3) k^' 



(2.6) 



(2.7) 



where for the convergence of the integral above is re- 
quested ub > — 3. We shall use to denote the magnetic 
field amplitude this quantity inst ead o f smearing the 
field at an additional scale A, as in iFPPi : 



{B')a = ^ / dkk^PB{k)e-'''~' 



(2.8) 



^ We choose the standard convention in which at present 
time toi ci{to) = 1- 

2 As Fourier transform and its inverse, we use: 
y(k,T) = / dxe^'^'^y (x, r) 



y(x,r) 



e-*"'y(k,T). 



where Y is a generic function. Note t hat we have changed 
our Fourier conventions with respect to lFPPj 
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The two definitions can however related simply by 



(uB + 3)r 



(2.9) 



The EMT of PMF is quadratic in the magnetic 
field and therefore its Fourier transfom is a convolu- 



tion. The two 
part of EMT ii 



{T:,(k)red(k')) = 



oint correlation function of the spatial 



dqdp 



647r5 

dqdp 

32^' 



In the appendices we explain the integration tech- 
nique and show the results for various spectral indexes 
n_B = 3,2,1,0,-1, —3/2 , —5/2: in this paper we add 
n_B = —5/2 to the previously studied ub > —3/2 values 
for scalar quantities studied in FPP. As for the scalar 
energy density and Lorentz force discussed in FPP, also 
the vector and tensor anisotropic stresses have support 
for Fourier modes with modulus smaller than 2kD- We 
show in Figure [1] the behaviour of scalar, vector and 
tensor quantities for ub ~ 2, —5/2, respectively. For 
B;(q)B;(k - q)B,„(p)B«|^^^-^)/2 the spectra for k <t: ko is: 



Ba (q) Bfc (k - q) (p)Bd (k' H^)^B^~ 



5/2 



We can then obtain scalar, vector and tensor correlation 
functions: 



whose slope, but not the a mplitude, agrees with 
iKahniashvih and Ratral l|2006l ). For riB = —5/2, 
|L(fc)|2 ~ (55/51)|p(fc)p for k <^ ko (for ub > -3/2 
we obtained \L{k)\^ ~ (ll/15)|p(fc)p for k <^ ko 



(2.14) 



^n*(S)(k)n(S)(k')) = <5a65cd{T:,(k)r,,(k')) 

(n*(^'(k)n(^)(k')) = fc„P.6(k)fc^P,d(k'){r:,(k)r,,(k')) 



(Ptc(k')Pw(k'; 



-PH(k')Pcd(k')){r:t,(k)red(k')): 



(v)r^^nW(k')) = fc„P,(k)fc^P,d(k'){r:,(k)r,,(k')) ""^^^ *™ contribution dom- 

■' inates over other ones in amplitude of the Fourier 

{n*j^'(k)n[f' (k')) = (Pia(k)Pj6(k) - ipj(k)P„6(k)) x spectra, in agreement with p revious numerical findings 

2 (|Brown and CrittendenlboOSi ). 

(2.10) The vector and tensor anisotropic stresses are 
shown with varying ub in Figure [S] The vector and 
tensor contributions have a fc-dependence very similar 
to the energy-density: for k <^ ko and ub > —3/2 
(ns = -3/2) |n(^)(/t)p , |n(^)(fc)|2 have a white noise 
(logarithmic divergent) spectrum with |n^"^'(fc)p ~ 
2|n'^'(fc)P, whereas for ub ~ —5/2 both become in- 
frared dominated as \U^'^\k)f ~ (94/25) IH^^' (it) p 



where Pi, 



-kikj. Such convolutions can be written 



in terms of spectra as follows: 



{n*'^'(k)m^'(k')) 



{n:'^'(k)n'.^)(k')) 



{n*j^)(k)n(p(k')) = 



(S) 



(fe)l'<5(k-k') 



n(^'(fc)|'p,(k)5(k-k' 



n'^'(fc)|^7Mijti(k)5(k — k') holds. The generic behaviour of |n^ iox k <^ ko 



and Ub > —3/2 is 



where Mijti = PitPji+ PuPjt - 
the spectra take the form: 



n(^'(fc)j^ 



2 7,2n+3 



A^k 



28 



(2.15) 



in 



PB{k)\^ 
(^'(fc)i^ 

n^^\k)f 



102471^ 
1 



dpPs(p)PB(ik-pi)(i(e.^i) 



5127r5 
1 



5127r5 



dpPi3(p)Pfl(|k-p|) > 

7')+7/3(m-7/3)] 

dpPi3(p)Ps(|k-p|) X 



(1 + 27^7^/3'), 



(2.12) 



(2.13) 



where 
P, 13 ^ 
agree. 



- p • (k - p)/lk - p|, 7 = k ■ 
k • (k — p)/lk — p|. These equations 
with in our Fourier convention, with p reviou s 
results by Mack. Kahniashvili and Kosowskvl l|2002l ). 
iDurrer. Ferreira and Kahniashvilil I 2000l ). One of the 
main results of our work is the calculation of the correct, 
i.e. without any approximation, expressions for these 
convolutions, given a power spectrum as in Eq.(2.5) 
with a sharp cut-off at ko- 

^ We use the convention that latin indexes run from 1 to 3 
while greek indexes run from to 3 



2567r4fc2"(3-h2ns) 15 ' 

The pole for ub = —3/2 in Eq. (|2.15p is replaced 
by a logarithmic diveregence in k in the exact result 
reported in the Appendix A; the resul t reported by 
iMack. Kahniashvili and Kosowskvl (|2002l ) has a factor 4 
instead of the factor 56/15 reported in Eq. (|2.15|l . Note 
also that the relation between the tensor and vector 
an isotropic stresses is different from t he on e reported 
in iMack. Kahniashvili and Kosowskvl (|2002i ). who pre- 
dict (in our conventions): \U<-'^'> {k)\'^ = \n'-^\k)\'^, 
this relation is obtained neglecting the angular part in 
Eas. (|2.12l2.13|l and is incorrect. 



3 PERTUBATIONS EVOLUTION WITH 
PMF 

The presence of PMF influences the cosmological per- 
turbations evolution mainly in three ways. PMFs carry 
energy and momentum at the perturbation level and 
therefore gravitate, influencing the metric perturba- 
tions. As a second point, they carry anisotropic stress, 
which adds to the ones already present in the plasma 
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0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 



Figure 1. Comparison of k'^\pg{k)\'^ (solid line), fc''|L(fc)p (large dashed line) fc^|n' (A:)|^ (small dashed line), fc''|n^^ C^)!^ 
(medium dashed line) in units of (iJ^}^/(10247r^) versus k/kjj. The left and right panel are for ub = 2 and ng = —5/2, 
respectively. 




given by neutrinos and photons, with the caveat that the 
photon anisotropic stress is negligible before decoupling 
epoch. Third, the presence of PMF induces a Lorentz 
force on baryons, which modifies their velocity. Due to 
the tight coupling between photons and baryons prior to 
the decoupling epoch the Lorentz force has an indirect 
effect also on photons during this regime. 

The evolution of metric perturbations is described 
by the Einstein equations. These are modified by the 
presence of PMFs that represent a source term as fol- 
lows: 

G^. = 87v{T^^ + T^^) , (3.16) 

where as usual t^^, represents the PMF EMT. The met- 
ric chosen in this work is: 

ds^ = a^(r) [-dr^ + [Si^ + hij)dx'dx^] (3.17) 

where hij can be decomposed into a trace part h and a 
traceless part consisting of its scalar, vector and tensor 



part (Ma and Bertschinger 1995): 

h^j = |<5,, + - ^V") M + 4 + 4 . (3.18) 

The vector part being constructed in terms of a diver- 
genceless vector hY 

hYj = dihj + djhY ■ (3.19) 
The tensor part is traceless and transverse {dihj^ = 0). 

4 THE SCALAR CONTRIBUTION 

We shall focus now on the magnetic scalar contribution 
to CMB anisotropics. The effect on metric perturbations 
is described by the Einstein equations with a source 
term given by the PMF EMT. We choose to work in the 
synchronous gauge where the scalar metric perturbation 
in the Fourier space is described by two scalar poten- 
tials, namely h{k, r) and rj{k, r). The Einstein equations 
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with the contribution of PMF in the synchronous gauge 
are: 

- = -47rGa^(E„ /9„5„ + pb) , 

h + 2TLh — 2k^Tj = — 87rGa^(E„ <?snpn5n 



+ 



5pE 



h + &ri + 2H{h + 677) - 2k^r) = -2'iTiGa x 
[E„(p„ + P„)a„ + crB], (4.20) 

where n represents the various species of the plasma, 
i.e. baryons, cold dark matter (CDM), photons and 
massless neutrinos. The conservation of the PMF 
EMT - V^r/;P^^ = - implies that Pi3(x,r) = 
Pb(x, ro)/a(r)* and the following relation between the 
magnetic anisotropic stress crs, the magnetic energy 
density pB and the Lorentz force L holds: 

aB = ^+L, (4.21) 

Such Lorentz force modifies the Euler equation for 
bar yons v elocity, leading to observational signatures 
fsee IFFpI for the most recent discussion on this effect). 



4.1 Initial Conditions 

The magnetized adiabatic mode initial conditions in the 
synchronous gauge deep in the radiation era are: 



36(15 + 4i?^) 



+ 



+ 



(-55 + 28P„)r2s 



336(15 + 4i?^) 1008(15 + 4i?^) 



I 4 4 
k T 



r\ — 2Ci — ^, ^ ^ ^ — ^ „ ^ Cifc r + 



6(15 + 4i?^; 
f2s(-55 + 28i?^) 



55Li 



168(15 + 4i?,.) 56(15 + 4i?^) 



, 2 2 
K r 



(5^ 

<5c 



-^B - lcik\^ + 
3 Ci 2 2 

'4 ~ T ^ ^ 

3 n C*! , 2 2 
--SZb —k T 

4 2 



+ 



6 ' 2(1 -Pi, 

»b(1-P.) ^ 
6P1, 



Lb 
2P, 



, 2 2 
K r 



SIb 3Lb 
~8" 8(1 -P.) 



Ci , 4 3 



+ 



SIb 3 Pb 
1~ 4(1-P, 

Pb 



k^T 



72 24(1 - P„) 



k\' 



+ 



(23 + 4P,) 
18(15 + 4P^) 

»b(1-P.) 
4P, 



+ 



3Lfl 
4 P, 



k^T 



(135 + 14P^)Lb 
84P^(15 + 4P„) 

(-270 + 161P^ + 2SRI)Qb 



504P^(15 



4Gi 



3(15 + 4P,) 

-!:!b(-55 + 28P^ 



4P.) 

Ab 
4P, 



+ 



56P^(15 + 4P^ 



+ 



3Ls 
' 4P, 

165Pb 
56P^(15 + 4P„) 



P3. = 



^B 3 Lb 
4P, 4 P, 



fcr , 



, 2 2 



(4.22) 



where Qb — Pb/{Pv + P7), Ls = L/(p^ + p^), 
Rv = pv/{pv + P7), Gi is the constant which charac- 
terizes the regular growing adiabatic mode as given in 
iMa and Bertschingerl (|l995l 'l. The resul t reported is dif- 
ferent from the one reported in IfppI because we are 
adding the first non-trivial terms in 0{Q,b ,Lb) for h , 5c 
and because of a different truncation in the neutrino hi- 
erarchy. The term C'(fc''r*) in /i (which also contains the 
next-to-leading term of the adiabatic mode) has been 
obtained by taking self-consistently the required order 
in all the variables in Eqs. (|4.22p . For simplicity we have 
written the 0(k'^T'^) term only in h (since it is the lead- 
ing term for the magnetic solution) and omitted these 
higher order terms in all the other variables. Here we 
also choose to truncate the h ierarchy at Fa,^ = in- 
stead of F'ii, = as in IFPPI This change affects the 
magnetic next to leading order terms in the velocity 
and anisotropic stress of neutrinos and in the metric 
perturbation rj. The equation for the evolution of the 
Ps^ is: 



F3. 



7 



koy . 



(4.23) 



while the neutrino anisotropic stress equation becomes: 

(4.24) 



Note how the presence of a SB of PMFs induces a 
new independent mode in matter and metric perturba- 
tions, i.e. the fully magnetic mode obtained by setting 
Gi = in Eq. (|4.22|) . This new independent mode is 
the particular solution of the inhomogeneous system of 
the Einstein-Botzmann differential equations: the SB of 
PMF treated as a stiff source acts indeed as a force term 
in the system of linear differential equations. Whereas 
the sum of the fully magnetic mode with the curva- 
ture one can be with any correlation as for an isocurva- 
ture mode, the nature of the fully magnetic mode - and 
therefore its effect - is different: the isocurvature modes 
are solutions of the homogeneous system (in which all 
the species have both background and perturbations). 
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whereas the fuUy magnetic one is the solution of the in- 
homogeneous system sourced by a fully inhomogeneous 
component. 

It is interesting to note how the magnetic contri- 
bution drops from the metric perturbation at leading 
order, although is larger than the adiabatic solution for 
photons, neutrinos and baryons. This is due to a com- 
pensation which nullifies the sum of the leading contri- 
butions (in the long-wavelength expansion) in the sin- 
gle species energy densities and therefore in the met- 
ric perturbations. A similar compensation exists for a 
network of topological defects, which does not carry a 
background energy-momentum tensor as the PMF SB 
studied her^B- 



4.2 Analytic Description of the Scalar 

Magnetic Contribution on Large Angular 
Scales 

In this subsection we give an analytic description for 
the scalar magnetic contribution to CMB anisotropics 
on large angular scales given by the initial conditions 
in Eq. (|4.22|l and computed by our modified version of 
CAME. The Sachs- Wolfe term is ST/T^sw ~ 5r/4 + ip, 
where 5r is the radiation density contrast and ij) is one 
of the metric potentials in the longitudinal gauge: 

ds^ = a^(?7) [-(1 + 2i))dTf + (1 - 2(t))'yijdx^dx^] (4.25) 

For adiabatic initial conditions is simply ST/T^sw ~ 
ip/3 since 5r — —8^p/3 in the matter dominated era on 
large scales. 

The Sachs- Wolfe term for the scalar magnetic mode 
generated by a SB of PMF can be obtained by using the 
initial conditions given in the synchronous gauge in Eq. 
(|4.22p . By making a gauge transformation, we obtain at 
leading order in the radiation era: 



Qb_ flBjl- R,) + 3Lb 
'4 15 + 4R^ 



14(15 + 4i..) ■ (^-26) 

In the radiation era, because of compensation, the met- 
ric potential are just proportional to Qb and not t o 
Q,B/k^ as obtained bv iKahniashvili and Rat"ral ()2006h . 
As we will show in the following, the same holds in the 
matter era. 

Now let us assume that in the matter era: 



5T _ Q.B 
T isw 4 



(4.27) 



and we compute the scalar contribution to CMB 
anisotropies for n > —3/2 by solving the integral: 

^ Note however that a network of topological defects does 
not scale with radiation and interacts only gravitationally 
with the rest of matter, i.e. a Lorentz term is absent. 



(4.28) 



where we have insterted an upper cut-off kn in order to 
use the infrared expansion of \pB{k)\, obtained in FPP: 



2i,2n-|-3 



\pB(k) 



A'k 



(4.29) 



1287r4fc2n(3 + 2ns) ' 
By using the result (|Abramowitz and Stegui]||l965l ) 

2 

2 y r 2 1 

dxxJe+i/2(x) = — [Je+i/2{y) - Ji-i/2{y)Ji+3/2iy)\ 



we obtain for n > —3/2 and y » 1: 

S ,magnetic ^ {n + 3)^(5^)^ 



(4.30) 



The scalar magnetic contribution to CMB anisotropies 
on large angular scales is therefore white noise 

^^S. magnetic ^ ^2) j^^. ^ ^ ^qq ^ _3/2. the 

slope in £ of Eq. (|4.30|l and of the numerical results ob- 
tained with our modified version of CAMB agree very 
well. By using an analogous procedure for ub = —5/2, 
we obtain analitically cf ''""'""'"'^ ~ 1/^ to be com- 



pared with the numerical results 



obtained 



with our modified version of CAMB. The a parameter in 
Eq. (|4.30|l we can fit from our numerical results inherits 
a dependence on ub , ko and is tipically ~ 0(0.01 — 0.1). 



5 THE VECTOR CONTRIBUTION 

In this section we shall review the evolution of vec- 
to r perturbatio ns induced by a SB of PMF as treated 
in iLewi j (|2004l ). The vector metric perturbation is de- 
scribed through a divergenceless vector: 



hYj = dihj -f djhi 



where we have 
dihi = 



(5.31) 



(5.32) 



The divergenceless condition assures that vector mode 
does not support density perturbations. The Einstein 
equations in the presence of PMF for the vector metric 
perturbations simply reduce to: 

+ 2Hh^ = -167rGa^(nl^' + H^^^ + ^T)/^ (5.33) 

Conservation equations for PMF also in the vector case 
reduce to a relation between isotropic and anisotropic 
pressure and the vector Lorentz force: 



v.p^ + v,n(p^ = Lf 



'j^Hj (5.34) 

The Lorentz force induced on baryons, in analogy with 
what we found for the scalar case, modifies the baryon 
vector velocity equation: 

Vb + Hvi, = ——{^neaGT{vh — «7) — — ) (5.35) 

Pb \0 p-y J 
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where we have neglected the baryon homogeneous pres- 
sure (pb/ Pb « !)• In order to investigate the effect of 
magnetized vector perturbations it is necessary to cal- 
culate the Fourier spectra for the vector projection of 
the PMF EMT and the Lorentz force. Since these two 
quantities are related by: 

= kU^P , (5.36) 

we need only one spectrum to compute for the vector 
part, as for the tensor part described in the next section 
and differently from the scalar part. 



and magnetic anisotropic stresses due to fact that mag- 
netic fields gravitate only at perturbative level. 

Note that the compensation between anisotropic 
stresses turns on only after neutrino decoupling, an 
epoch which is much earlier than the usual initial time 
at which cosmological perturbations are evolved in an 
Einstein-Boltzmann code. The detailed study of the evo- 
lution of gravitational waves during neutrino decoupling 
is an interesting issue, but clearly beyond the purpose 
of the present project. 



6 THE TENSOR CONTRIBUTION 

Inflationary tensor modes, namely primordial gravita- 
tional waves, are a key prediction of the standard infla- 
tionary model, and therefore their indirect observation 
through CMB anisotropies is one of the crucial point of 
modern cosmology. However PMF carrying anisotropic 
stress are themselves a source of tensor perturbations. 
Therefore the presence of PMF affects inflationary ten- 
sor modes and creates a new independent fully magnetic 
tensor mode in analogy to what we found for scalar per- 
turbations. The evolution equation for the metric tensor 
perturbation hij is: 



hi 



+ 2Uh,, -f k'h^, = 167rGa"(p,<, + ') ■ (6.37) 

which, for each polarization state deep in radiation era 
reads 



'i.T , ^ ; T 
hk + -hk ■ 



where 11^ represents the time independent variable 
n^'/pT. The large scales solution to this equation can 
be found expanding hh in powers of (fcr). In order to 
keep the leading and the next-to-leading terms we need 
to take into account the neutrino octopole (J3), trun- 
cating the propagation of anisotropic stress through 
higher moments by posing J4 = 0. Hence the neutrino 
anisotropic stress evolves according to 



(T) 



i?^)n 



fT(T)l 



(6.38) 



■ST) 



--hu - -J3 



and the solution is then: 

5(fcr)2 



hk 



JT) 



A 



1 



(1 



2(15 -I- 4J?, 



- ■R^) j^(T) 



+ 



15(1 



(6.39) 



R^)Il^P{krf 



14(15 
15(fcr)2 



4i?,) 



14(15 + 4i?,) 



+A 



2(fcr)2 



3(15 -I- 4i?,) 



(6.40) 



The presence of magnetic fields is responsible for 



the new leading term in al - otherwise absent. This 
is the so-called compensation between coUisionless fiuid 



7 RESULTS FOR CMB ANISOTROPIES 

In this section we now present the temperature and po- 
larization CMB spectra including the full contribution 
of SB of PMF. In addition to the Ce obtained by the adi- 
abatic mode in absence of primordial magnetic fields, we 
add the three contributions C^'^'^ described in Sects. 
4,5,6, i.e. scalar, vector, tensor, respectively, computed 
separately by our modified version of CAMB. In Fig. 3 
and Fig. 4 we show the results for ub = 2 and —5/2, re- 
spectively. For the initial conditions of the scalar mag- 
netic mode, we use as initial conditions Eq. (4.22) with 
Ci — 0; for those of the tensor mode the ones in Eq. 
(6.40) with A — 0. For the initial conditions for the 
vector mode we use the ones alr eady implemented in 
CAMB, described in lLewisI (|2004l ). AU the formulae for 
magnetic spectra needed - Ub , Lb , Hv, Ht - are given 
in our appendices; the si gns of Lb and Sis are taken as 
opposite, as explained in FPPl . 

For all the values of ns considered here, the CMB 
temperature pattern generated by the SB of PMF is 
dominated by the scalar contribution at low and inter- 
mediate multipoles; the vector contribution takes over 
the scalar one at high multipoles, whereas the tensor one 
is always subleading with respect to scalar and vector. 

It is interesting to note that the B polarization sig- 
nal due to the vector contribution is always larger than 
the tensor one. The B mode produced by vector per- 
turbations has a power spectrum which can be steeper 
than the one produced by lensing with a peak around 
£ ~ few X 10^; therefore, for suitable values of the mag- 
netic field amplitude the B mode produced by a SB of 
PMF can be larger than the lensing one for any ub- Fig. 
5 shows how the vector contribution to the B spectrum 
depends on ub- For ub > —3/2 the B spectra from the 
vector contribution are almost indistinguishable for dif- 
ferent Ub, because H^' is white noise for k ko', for 
—3 < Ub < —3/2 we note a dependence of the B spec- 
trum on ns. Analogous dependence on ub also holds 
for the vector contribution to TT. 
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10 100 'UUO 10 100 1000 
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Figure 3. CMB anisotropies angular power spectrum for temperature (TT hereafter, top-left panel), temperature-E polar- 
ization cross correlation (TE hereafter, top-right panel), E polarization (EE hereafter, bottom-left panel), B polarization (BB 
hereafter, bottom-right panel). The solid line is the adiabatic scalar contribution in TT, TE, EE panels, whereas it is the tensor 
homogeneous contribution in the BB panel (for a tensor-to-scalar ratio r = 0.1); the dotted, dot-dashed, dashed are the scalar, 
vector and tensor contributions of a SB of PMF respectively for ^J(B^ = 7.5 nG, ng = 2 and ko = 27rMpc~^. The dotted 
line in the BB panel is the lensing contribution. The cosmological parameters of the flat KCDM model are Cj, = 0.022, 
n^h'^ = 0.123, Zro = 12, Us = 1, Ho = 100/ikms~i Mpc"! = 72kms-l Mpc"!. 



8 CONCLUSION 

We have obtained the Fourier spectra of the rele- 
vant scalar, vector and tensor energy-momentum com- 
ponents of th e SB of PMF extending the method 
used in IFPPI only for the sc alar. As already dis- 
cussed for the scalar sector in IFPPI . we have shown 
how the correct evaluation of the convolution inte- 
grals leads to differences in the vector and tensor 
p arts of the PMF Fourier spectrum previo usly found 
in lMack. Kahniashvili and Kosowskvl (|2002l ). 

We have then shown the comparison of the scalar, 
vector and tensor contributions to CMB anisotropies of 
the new inhomogeneous modes generated by the SB of 
PMF, by using the correct convolutions for the energy- 
momentum tensor of the PMF SB. We have shown that 
the dominant contributions are from scalar and vector 
perturbations, respectively for low and high £. We have 



given an analytic description of the Sachs- Wolfe contri- 
bution of the scalar mode, which agrees very well with 
the numerical result of our modified version of CAMB 
(|Lewis.Challinor and Lasenblvll2000l ) and takes into ac- 
count the compensation effect on large scales which is 
generic when a fully inhomogeneous source is present. 
The slope in £ of the vector power spectrum we ob- 
tain num ericall y agrees very well with pre vious analytic 
(iMack. K ahnias hvili and Kosowskvl |2002| ) and n umeri- 
cal ( LewisI |2004) results. As already found in l|Lewij 
12004 ). the B signal by vector perturbations, with a 
peak slightly dependent on ns around £ ~ 2000, has 
a shape difi'erent either from the infiationary gravita- 
tional waves or the lensing signal. We have character- 
ized its dependence on ns and shown how this is non 
trivial f or hb < —3/2: such signal can be constrained by 
Planck ( | [Planck Collaboration]! |2006|) and future small 
scale CMB polarization experiments. 
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Figure 4. CMB angular jjowcr sijcctrum for TT (left top panel), TE (left top panel), EE (bottom left), BB (bottom right). The 
solid line is the adiabatic scalar contribution in TT, TE, EE panels, whereas it is the tensor homogeneous contribution in the 
BB panel (for a tensor-to-scalaj: ratio r = 0.1); the dotted, dot-dashed, dashed are the scalar, vector and tensor contributions 
of a SB of PMF respectively for y^{B2) = 7.5 nG, ub = —5/2 and kr> = 27rMpc ^. The dotted line in the BB panel is the 
lensing contribution. The cosmological parameters of the flat ACDM model are the same as in Fig. 3. 
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10 100 1000 



Figure 5. Vector contributions to the CMB angular power spectrum for BB. The soUd line is the tensor homogeneous contri- 
bution for a tciisor-to-scalar ratio r = 0.1 and the dotted line is the lensing contribution with cosmological parameters as in the 
previous figures; the triple-dotted, long dashed, dot-dashed and dashed are the vector spectra obtained with (B^) = 7.5 nG, 
ku = 27rMpc~^, for ub = —2.5, —1.5, —1, 2, respectively. Note how the spectra for ng = —1 and ub = 2 are super-imposed 
since the Fourier spectra of the vector part of the PMF EMT are both white noise for fe <S feo for these spectral indexes. 
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APPENDIX A: EMT FOURIER SPECTRA 

We use the convolutions for the PMF EMT spectra with the parametrization for the magnetic field power spectrum 
given in Eq. (|2.5p . Since Ps(fc) = for k > ko, two conditions need to be taken into account; 

p < ko , |k - p| < fcu . (Al) 

The second condition introduces a fc-dependence on the angular integration domain and the two allow the energy 
power spectrum to be non zero only for < A; < 2fc_D. Such conditions split the double integral (over 7 and over p) in 
three parts depending on the 7 and p lower and upper limit of integration. For simplicity we normalize the Fourier 
wavenumber to ko and we show the integrals with this convention. A sketch of the integration is thus the following: 

1) < fc < 1 

1 /.I pi i-l-k nl 

...+ I dp I ^ ^ dj ...= I dpla{p, k) + I dph{p, k) 

— 4f 



1 Jl-k ~ Jo Jl-k 



2) 1 < fc < 2 

dp [ d-f ...= I dph[p,k) (A2) 



Particular care must be used in the radial integrals. In particular, the presence of the term \k — p]"'*'^ in both 
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integrands, needs a further splitting of the integral domain for odd n: 

with p < k 

fc 

(i-fe) 



(i-fc) 



dp - 



dp 



(i-fc) 



I dp... with p>fc 



fc>l/2 
fc<l/2 



with p < k 

with p > k 
with p < k 
with p > k 



I dp^ I l<k<2 dp... with p < fe 



Following the scheme (|A2p we can now perform the integration over p. Our exact results are given for particular 
values of n_B. 
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APPENDIX B: SCALAR PART OF THE LORENTZ FORCE 

In order to compute the scalar contribution of a SB of PMFs to the cosmological perturbations, the convolution for 
the scalar part of the Lorentz Force power spectrum is also necessary. The scalar anisotropic stress can be obtained 
directly from its relation with the Lorentz force and the magnetic energy density in Eq. (|4.2H) . We report here the 
result for the Lorentz force convolution: 

\mf - io2L^a^ J '^'^ PB(p)Pfl(|k-p|)[l + M2+ 47/3(7/3 (B1) 
with the magnetic field power spectrum in Eq. p.5|l for particular values of n_g. 



Bl riB^S 



44 _ 2fc I 556P _ 4p_ , 164fc'* , 4k'' _ llfc^ r () < h < ^ 

135 3 ^ 735 9 ~'~ 1575 ~'~ 2079 11025 U fi, 1 

for 1 < fc < 2 



421.9 135 3 735 9 1575 207' 

\r(U]\^ _ ) 44 , 64 16 I 88 1 2k 

1 V ;inB=3 5127r4fc6 S 135 24255fe5 ,945fe3 T„ 525fe 3 

556*;'' I 4p _ 164*" _ Ak" , 11*:^ 
^ 735 9 1575 2079 33075 



B2 



riB 



\m\ls=2 



A^kjy 


' 44 


2k 8k^ 




P 13fc'' 


5127r4fc4 


105 


3 ' 15 


6 


240 ' 6720 



© 0000 RAS, MNRAS 000, 000-000 



The contribution of PMFs to CMB anisotropies 15 



B3 



riB 



l^(fc)in = l = 



44 _ 2fc I 32fc-' I 4k2_ _ 
75 3 105 315 25 

_44 I 64 16 J 8^ I 2fc 

75 + 1575*!= 105*3 ISfe 3 



32fc^ 
105 



for < fc < 1 
+ for l<fc<2 



B4 ns=0 



5127r4 
log 11 



43 
48 



1 



+ 



16fc4 32fe3 
log j 1 - fc 
8fc 



+ 



+ ■ 



+ 



48A;2 
P log|l 



13 
192fc 



67fc P 
"96" ^ 48 



17 P 

384 



log|l 



16fe5 



B5 = -1 



|-^(^)lni5=-l 



5127r4 



15 105 



44 
" 15 



64 I 16 
105*5 15*3 



I 8 I 2fc I 4k± 
3fc 3 ^105 



for < fe < 1 
for 1 < fc < 2 



B6 



riB 



-3/2 



\m\l- 



-3/2 



5127r4 



1755V 1-* 
I 88 _ 
15* 



2925*5 



2925V l-**5 
32* 



222L 4- 128 
15 135fe3 



1775V1-** 1755\ 
10616 2048 



1-* 

2048 



1775V 1-* 



+ 



8775 V l-fc fc3 
log[l+\/l-fc] 44 log fc 



1755\/-l-l-fe 2925fe5 2925\/-l-|-fefe5 2925\/^l+I*'* 
I 128 I 9088 3776 L 88 , 10136 



135*3 
32* 



8775\/-l+**3 8775\A^1+Ifc' 
44 arctan 



1755\/-l + fefe 



4- 



64fe^ 



1775V-1-H* 1775V-l-Hfe 



H-fe 



for < fc < 1 



for 1 < fc < 2 



B7 UB = -5/2 

l^(fc)lnB = -5/2 = 



A^kl 



5127r4fc|, 



32 



+ 



2048 



2048 



1024 



1155vjl-fcl 1155fc5 1155\/|1 - k\P 1155\/|1 - ^1^" 
1664 12976 88 13648 64fc 



+ 



1155v/ll - fc|fc^ 1155^/11 - k\P 



+ 



15fc 



— + ■ 



128 
105fc3 



1155^/11 - fcjfc 1155v/|l - fc| 



Acknowledgements. We wish to thank Chiara Caprini, Ruth Durrer, Antony Lewis and Kandaswamy Subra- 
manian for comments and discussions. This work has been done in the framework of the Planck LFI activities and is 
partially supported by ASI contract Planck LFI Activity of Phase E2. 



REFERENCES 

Abramowitz M. and Stegun I., Handbook of mathematical functions with formulas, graphs, and mathematical table, 

New York: Dover Publishing (1965). 
Brown L and Crittenden R., Phys. Rev. D 72 (2005) 063002 

Caprini C, Durrer R. and Kahniashvih T., Phys. Rev. D 69 (2004) 063006 

Durrer R., "Cosmic Magnetic Fields and the CMB", New Astron. Rev. 51 (2007) 275 [arXiv:astro-ph/0609216] . 

Durrer R.,Ferreira P.G. and Kahniashvih T.,Phys. Rev. D 61 (2000) 043001 

Finelh F., Pad F., Paoletti D., 2008, Phys. Rev. D, 78,(2008) 023510 

Giovannini M. and Kunze K.E.,Phys. Rev. D77 (2008) 063003 

Kahniashvih T. and Ratra B., Phys. Rev. D 75,(2006) 023002 

Koh S., and Lee C.H., Phys. Rev. D 62, (2000) 083509. 



© 0000 RAS, MNRAS 000, 000-000 



16 D. Paoletti, F.Finelli,F.Paci 



Lewis A., Phys. Rev. D, 70,(2004) 043011 

Lewis A., Challinor A. and Lasenby A., Astrophys. J. 538 (2000) 473 
Ma C. P. and Bertschinger E., Ap. J., 455 (1995) 7. 

Mack A., Kahniashvili T., Kosowsky A., Phys. Rev. D, 65,(2002) 123004 

[Planck Collaboration], "Planck: The scientific programme," arXiv:astr o-ph/0604069| (2006). 
Seshadri T. R and Subramanian K., Phys. Rev. Lett. 87 (2001) 101301 
Subramanian K., "Magnetizing the Universe", arXiv:astro-ph/0802.2804 [astro-ph] 
Subramanian K. and Barrow J. D., Phys. Rev. Lett. 81 (1998) 3575 
Yamazaki D. G., Ichiki K., Kajino T., and Mathews G.J., apj, 646, 719 (2006). 



© 0000 RAS, MNRAS 000, 000-000 



